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Pionic couplings g and g in the static heavy quark limit 



1. Physics Motivation 



Static quark limit of QCD offers a simplified framework to solving the non-perturbative dy- 
namics of light degrees of freedom, which is highly important for various phenomenological studies 
of weak interactions involving the heavy-light mesons. In the exact niQ — » oo limit, the heavy quark 
symmetry (HQS) constrains the QCD dynamics of light quarks and soft gluons in the heavy-light 
hadrons to be invariant under the change of the heavy quark flavour and/or its spin Uq). As a result 
the total angular momentum of the light degrees of freedom becomes a good quantum number (n), 
and therefore the physical heavy-light mesons come in mass-degenerate doublets: 



[D(0-),D*(1-)}, [D* Q (0 + ),D[(n], [D 1 (1 + ),D^(2+)] 
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where on the example of charmed states we remind the reader of spectroscopic labels. HQS is a 
basis for the systematic expansion of the QCD Green functions in power series of Aqcd/'mq, known 
as the heavy quark effective theory (HQET). At each order, one has to solve the non-perturbative 
QCD dynamics of the light degrees of freedom which can actually be done only by the numerical 
QCD simulations on the lattice. 

Due to inability to simulate directly the dynamics 
of light quarks that are close enough in mass to the 
physical u- or J-quark, one has to make various chi- 
ral extrapolations, which induce large systematic errors 
to the lattice results (a typical example is the one of 
SU (3) flavour breaking ratio of the pseudoscalar me- 
son decay constants, /bJ/bj)- In addition the insuffi- 
cient size of the lattice box with the periodic bound- 
ary conditions afflicts the propagation of the light pi- 
ons on the lattice. To have a good handle on these 
two problems one can rely on the effective theory that 
describes the dynamics of light degrees of freedom in 
terms of pseudo-Goldstone bosons (pions,, for short), 
usually referred to as heavy meson chiral perturbation 
theory (HMChPT) HMChPT is suitable when work- 
ing with (nearly) massless quarks and thus it is comple- 
mentary to the light quark dynamics that is directly ac- 
cessible from the lattice. Similar to the standard ChPT, 
where the pion axial coupling (f n ) is a parameter of the 
theory, in the HQChPT Lagrangians the axial couplings 
of a pion to doublets of heavy-light mesons are param- 
eters of the theory. In other words, they have to be fixed 
elsewhere. The most well-known such a coupling is the 
one with the (l/2)~ -doublet, known as g, whose value 

was experimentally established in the case of heavy charm quark 0. That value appeared to be 
much larger than the ones predicted by most of the relativistic quark models and by all the QCD 
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Figure 1: Results for g with the heavy 
quark being static, bottom or charm. As- 
terisk on the lattice result for g/, indicates 
that the result is not obtained after interpo- 
lating between the static and charm heavy 
quark cases. The numbers are taken from 
refs. b| Q g], in descending order. 
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sum rules. In fig. [j], we show the panorama of the results for that coupling obtained by using lat- 
tices and light cone QCD sum rules (LCSR). Lattice values are obtained both in the static heavy 
quark limit []|, ^] and with the heavy quark of the mass around the physical charm quark mass |^]. 
However, all the lattice computations of that coupling are obtained in the quenched approximation. 
In this note, we present the first unquenched result for g, with Nf = 2 dynamical light quarks and 
in the static heavy quark limit (indicated as "New" in fig. [|. In addition, we report the first result 
for g, the coupling of the soft pion to the (l/2) + -doublet of heavy-light mesons. 



2. Definitions and Correlation functions to be computed 

In the limit in which the heavy quark is infinitely heavy and the light quarks are massless, the 
axial couplings of the charged pion to the heavy-light mesons, g and g, are denned via 

(H\J\H*(e))=ge, (H^\H[(e))=ge, (2.1) 

where the non-relativistic normalisation of states \Hq*^) is assumed, (Hi(v)\Hj(v')} = 8jj8(v — v'). 
All heavy-light hadrons are at rest (v = v 1 = 0), and therefore the soft pion which couples to the 
axial current, ^ = uy^d, is also at rest, q = 0. is the polarisation of the vector and the axial 
heavy-light meson respectively. 

The standard strategy to compute the above matrix elements consists in evaluating the follow- 
ing two- and three- point conelation functions: 
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(2.2) 



where (...) v denotes the average over independent gauge field configurations. The gauge field 
configurations used in this work are the ones produced by the SPQcdR collaboration [j7p. The 
interpolating fields that we use to produce the desired (1/2) T heavy-light mesons are: P = hy^q, 
Vi = hjiq, S = hq, and Aj = hjij^q, with q being either u- or J-quark field. In eq. ( |2.2| ) we also 
express the correlation functions in terms of quark propagators: light ones, y u j(0,x), and the 
infinitely heavy one (the Wilson line), 



T=0 



(2.3) 



where Uo is the temporal component of the link variable. The spectral decomposition of the three 
point function Cs(t x ) reads 

C 3 (t x ) = £ [^0]P|^e-^<fl»|^|ff;)e-^- £ *>*^|J2l^] , (2.4) 

Z„ Z,„£j(m) 
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with the sum running over the ground states (m = n = 0) and the radial excitations (m, n>0) which 
are heavier and thus exponentially suppressed. If only the diagonal terms were different from zero 
[i.e., the terms in the sum (2.4) with m = n], then the function C 3 (t x ) would be ^-independent 
because the HQS ensures that En* — En . This is indeed what one observes from the lattice data: 
even for very small values of t x the correlation functions C 3 (t x ) and C$(t x ) exhibit plateaus (c.f. 
fig. |J). However, even after leaving out the non-diagonal terms (m ^ n), we still have a problem 
of contamination by the radial excitations (m = n > 0) in eq. (j2"!4|). The latter are taken care of by 
implementing the smearing procedure which highly enhances the overlap of the interpolating fields 
with the lowest energy states. Once that is been done, it is then trivial to extract the couplings g 
and g, namely 



C 3 (t x ) 



•(ml), 



C 3 (t x ) 
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(2.5) 



where (Zq" 1 ) 2 is extracted from the fit to C2{t) = ^|Z^ m ' | 2 exp (—E n t), with smeared sources. Like- 

n 

wise (Z™ ) 2 is extracted from C^/). Before going to the results, we should stress again that the 
above identification with the couplings g and g is valid only in the chiral limit, g = lim g(m 2 n ). 
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Figure 2: Correlation functions C3 (t x ) and C3 (t x ) in 
which the sources are fixed at and t y = 12. 



3. Lattice details and results 

Our dataset consists of 50 gauge field con- 
figurations which include Nf = 2 dynamical quarks 
whose masses, with respect to the physical strange 
quark mass, are within the range 0.7 < m q /m% hys < 
1.3. We use the (unimproved) Wilson quarks 
and work on 24 3 x 48 lattices at j3 =5.8, which 
corresponds to the lattice spacing of a « 0.06 fm, 
and the physical volume (1.5 fm) . The cor- 
responding "sea" quark hopping parameters (fc sea ) 
are listed in table [[]. More information about 
the simulations can be found in ref. [pp. In all 
our computations we keep the valence and sea 
light quarks equal. 

In the computation we used the Eichten- 
Hill static quark action [§J, supplemented by 
^] which is known to substantially improve the 



the hypercubic blocking procedure (HYP) of ref. 
signal/noise ratio in the correlation functions (|2~2|). In that way we were able to monitor the the 
efficiency of the smearing, i.e., by confronting the signal of Czit) by using both local sources with 
the signal in which both sources are being smeared. The generic local source operator is smeared 
as in ref. [ffO], h(x)Yq(x) — > h(x)Yq s (x) where, 
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U F being the so called fuzzed link variable. The parameters R = 3 and R max = 5 are tuned to 
enhance the overlap with the lowest lying state, which we were able to check from the analysis of 
the two point functions C 2 (t) and C 2 (t). We obtain that Z* m /Z l ° c - » 150, while Z\ m /Z l ° c - < 0.05. 
Similar situation holds for the scalar-scalar correlation function, and it is therefore safe to say that 
the pollution of our results for g{m 2 n ) and that comes from the transitions among higher 

excitations is to a huge extent suppressed by the smearing procedure. 
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Table 1: Results for the orbital splitting Eq — Eq, the couplings g and g, from the unquenched Nf = 2 
simulations in which the dynamical quark is of the mass corresponding to the pion whose masses are also 
given. 



Our results for g{m\) and g{m\) are given in table [j], and plotted in fig. [| These results are 
obtained by fixing one of the source operators to t y = 12. The stability of these values is checked by 
using t y = 11, and t y = 13. In the same table we also give the mass difference between the lowest 
lying states belonging to (l/2) + and (l/2)~ doublets. On the basis of our results, we make the 
following observations: 
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Figure 3: Our new (unquenched) results against the old (quenched) ones. The dashed lines correspond to 
the range of values that g(0) and ^(0) can assume after the chiral extrapolation. 
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• The values of g(m„) from the unquenched simulations are consistent with the quenched 
values we reported in ref. [||] (see also fig. |||); 

• The values of g{m\) are clearly negative and their absolute values much smaller than those 

ofg(m|); 

• The statistical and systematic quality of our data is not good enough to discuss the subtleties 
of chiral extrapolations. Instead, we quote 0.4 < g < 0.6, and —0.3 < g < —0.1, as our 
current estimates; 

• Less prone to statistical noise is the ratio g/g, for which we obtain 

g/g = -0.34(3), -0.43(9), -0.39(7), -0.31(6), (3.2) 
in the order of decreasing light quark mass; 

• Our results are in a clear conflict with the assumption, g = g, that stems from the parity 



doubling models [11]; 



• Our values of the orbital splitting E$ — Eq, together with the couplings g and g do not help 
to explain the puzzling experimental phenomenon that {Eq — Eq) s < (Eo — Eo) u ,d [|2|]- To 
address this question more seriously it is essential to reduce the sea quark masses, work on 
larger lattice volumes and increase the statistical quality of our data. 

References 

[1] R. Casalbuoni, A. Deandrea, N. Di Bartolomeo, R. Gatto, F. Feruglio and G. Nardulli, Phys. Rept. 
281 (1997) 145 [hep-ph/9605342]. 

[2] A. Anastassov et al. [CLEO Collaboration], Phys. Rev. D 65 (2002) 032003 [hep-ex/0108043]. 

[3] A. Abada et al, Phys. Rev. D 66 (2002) 074504 [hep-ph/0206237]; Nucl. Phys. Proc. Suppl. 119 
(2003)641 [hepdat/0209092]. 

[4] A. Khodjamirian et al, Phys. Lett. B 457 (1999) 245 [hep-ph/9903421], see also P. Ball and 
R. Zwicky Phys. Rev. D 71 (2005) 014015 [hep-ph/0406232]. 

[5] A. Abada et al, JHEP 0402 (2004) 016 [hepdat/03 10050]. 

[6] G. M. de Divitiis et al [UKQCD Collaboration], JHEP 9810 (1998) 010 [hepdat/9807032]. 
[7] D. Becirevic et al, hepdat/0509091, and hepdat/05 10014. 
[8] E. Eichten and B. Hill, Phys. Lett. B 234 (1990) 511. 

[9] A. Hasenfratz and F. Knechtli, Phys. Rev. D 64 (2001) 034504 [hepdat/0 103029]. 

[10] P. Boyle [UKQCD Collaboration], J. Comput. Phys. 179 (2002) 349 [hepdat/9903033]. 

[11] W. A. Bardeen, E. J. Eichten and C. T. Hill, Phys. Rev. D 68 (2003) 054024 [hep-ph/0305049]; 
M. A. Nowak, M. Rho and I. Zahed, Acta Phys. Polon. B 35 (2004) 2377 [hep-ph/0307102]. 

[12] T. Mehen and R. P. Springer, Phys. Rev. D 72 (2005) 034006 [hep-ph/0503134]; D. Becirevic, 
S. Fajfer and S. Prelovsek, Phys. Lett. B 599 (2004) 55 [hep-ph/0406296]. 



212/6 



